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Nomenclature 


2D: 

2 dimensional 

4D: 

4 dimensional 

ADS-B: 

Automatic Dependence Surveillance Broadcast 

CAS: 

Calibrated Airspeed 

DTG: 

Distance-To-Go 

MSL: 

Mean Sea Level 

STAR: 

Standard Terminal Arrivals 

TAS: 

True Airspeed 

TCP: 

Trajectory Change Point 

TOD: 

Top-Of-Descent 

TTG: 

Time-To-Go 

VTCP: 

Vertical Trajectory Change Point 

Subscripts 


Subscripts associated with waypoints and TCPs, e.g., TCP 2 , denote the location of the waypoint or TCP in 
the TCP list. Larger numbers denote locations closer to the end of the list, with the end of the list being 
the runway threshold. Subscripts in variables indicate that the variable is associated with the TCP with 
that subscript, e.g., Altitude 2 is the altitude value associated with TCP 2 . 

Units and Dimensions 

Unless specifically defined otherwise, units (dimensions) are as follows: 


time: 

seconds 

position: 

degrees, + north and + east 

altitude: 

feet, above MSL 

distance: 

nautical miles 

speed: 

knots 

track: 

degrees, true, beginning at north, positive clockwise 
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Abstract 


This document describes an algorithm for the generation of a 
four dimensional trajectory. Input data for this algorithm are 
similar to an augmented Standard Terminal Arrival (STAR) with 
the augmentation in the form of altitude or speed crossing 
restrictions at waypoints on the route. This version of the 
algorithm accommodates descent Mach values that are different 
from the cruise Mach values. The algorithm calculates the 
altitude, speed, along path distance, and along path time for 
each waypoint. Wind data at each of these waypoints are also 
used for the calculation of ground speed and turn radius. 


Introduction 

Concepts for self-spacing of aircraft operating into airport terminal areas have been under 
development since the 1970's (refs. 1-20). Interest in these concepts has recently been renewed 
due to a combination of emerging, enabling technology (Automatic Dependent Surveillance 
Broadcast data link, ADS-B) and the continued growth in air traffic with the ever increasing 
demand on airport (and runway) throughput. Terminal area self-spacing has the potential to 
provide an increase in runway capacity through an increase in the accuracy of runway threshold 
crossing times, which can lead to a decrease of the variability of the runway threshold crossing 
times. Current concepts use a trajectory based technique that allows for the extension of self- 
spacing capabilities beyond the terminal area to a point prior to the top of the en route descent. 

The overall NASA Langley concept for a trajectory-based solution for en route and terminal 
area self-spacing is fairly simple and was originally documented in reference 21. By assuming a 
4D trajectory for an aircraft and knowing that aircraft’s position, it is possible to determine where 
that aircraft is on its trajectory. Knowing the position on the trajectory, the aircraft’s estimated 
time-to-go (TTG) to a point, in this case the runway threshold, is known. To apply this to a self- 
spacing concept, a TTG is calculated for a leading aircraft and for the ownship. Note that the 
trajectories do not need to be the same. The nominal spacing time and spacing error can then be 
computed as: 

nominal spacing time = planned spacing time interval + traffic TTG. 

spacing error = ownship TTG - nominal spacing time. 

The foundation of this spacing concept is the ability to generate a 4D trajectory. The algorithm 
presented in this paper uses as input a simple, augmented 2D path definition (i.e., a traditional 
STAR, with relevant speed and altitude crossing constraints) along with a forecast wind speed 
profile for each waypoint. The algorithm then computes a full 4D trajectory defined by a series of 
trajectory change points (TCPs). The input speed (Mach or CAS) or altitude crossing constraint 
includes the deceleration rate or vertical angle value required to meet the constraint. The TCPs 
are computed such that speed values, Mach or CAS, and altitudes change linearly between them. 
TCPs also define the beginning and ending segments of turns, with the midpoint defined as a fly- 
by waypoint. The algorithm also uses the waypoint forecast wind speed profile in a linear 
interpolation to calculate the wind speed at the altitude the computed trajectory crosses the 
waypoint. Wind speed values are then used to calculate the ground speeds along the path. 

The major complexity in computing a 4D trajectory involves the interrelationship of ground 
speed with the path distance around turns. In a turn, the length of the estimated ground path and 
the associated turn radius will interact with the waypoint winds and with any change in the 



specified speed during the turn, i.e., a speed crossing-restriction at the waypoint. Either of these 
conditions will cause a change in the estimated turn radius. The change in the turn radius will 
affect the length of the ground path which can then interact with the distance to the deceleration 
point, which thereby affects the turn radius calculation. To accommodate these interactions, the 
algorithm uses a multi-pass technique in generating the 4D path, with the ground path estimation 
from the previous calculation used as the starting condition for the current calculation. 


Algorithm Overview 

The basic functions for this trajectory algorithm are shown in figure 1. Note that waypoints are 
considered to be TCPs but not all TCPs are waypoints. 

For the 2D input, the first and last waypoints must be fully constrained, i.e., have both a speed 
and altitude constraint defined. With the exception of the first waypoint, which is the waypoint 
farthest from the runway threshold, constraints must also include a variable that defines the 
means for meeting that constraint. For altitude constraints, this is the inertial descent angle; for 
speed constraints, it is the air mass CAS deceleration rate. A separate, single Mach / CAS 
transition speed (CAS) value may also be input for profiles that involve a constant Mach / CAS 
descent segment. Additionally, an altitude / CAS restriction (e.g., in the U.S., the 10,000 ft / 250 
kt restriction) may also be entered. 

The algorithm computes the altitude and speed for each waypoint. It also calculates every point 
along the path where an altitude or speed transition occurs. These points are considered vertical 
TCPs (VTCPs). TCPs also define the beginning and ending segments of turns, with the midpoint 
defined as a fly-by waypoint. Turn data are generated by dividing the turn into two parts (from 
the beginning of the turn to the midpoint and from the midpoint to the end of the turn) to provide 
better ground speed (and resulting turn radius) data relative to a single segment estimation. A 
fixed, average bank angle value is used in the turn radius calculation. The algorithm also uses the 
forecast wind speed profile for a waypoint in a linear interpolation to calculate the wind speed at 
the altitude the computed trajectory crosses the waypoint (if the crossing altitude is not at a 
forecast altitude). For non-waypoint TCPs, the generator uses the forecast wind speed profile 
from the two waypoints on either side of the TCP in a double linear interpolation based on 
altitude and distance (to each waypoint). Of significant importance for the use of the data 
generated by this algorithm is that altitude and speeds (Mach or CAS) change linearly between 
the TCPs, thus allowing later calculations of DTG or TTG for any point on the path to be easily 
performed. 
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Initialize the waypoint turn data : Waypoints that have more than a 3 degree change in 
ground track from the previous waypoint are considered turn-waypoints. Mark each as 
a turn-waypoint and insert a turn-entry and turn-exit TCP on each side of this waypoint. 


while this is not the last loop 


Compute the TCP altitudes : Beginning at the runway (the last waypoint) work 
backwards and compute the altitude at each prior TCP. If an altitude is computed to be 
reached prior to the previous TCP, insert a new altitude TCP. 


Copy crossing angles : Beginning at the runway, for TCPs that do not have crossing 
angles, copy the downstream angle into this TCP. 


Do Descent Mach: If the descent Mach is different than the cruise Mach, change the 
descent waypoint Mach as necessary and add any required, additional TCPs. 


Compute the Mach / CAS TCP: If required, compute the Mach / CAS transition 
altitude. Compute the DTG to this altitude and insert the Mach / CAS TCP. 


if flag for an altitude / CAS 
restriction, based on prior 
computations, is 
\ true 


Compute an altitude / CAS TCP : In the U.S., this would be the 1 0,000ft / 250kt speed 
restriction. If the speed crossing the trajectory at the specific altitude is greater than the 
CAS restriction (from the test on the first time through), place a speed restriction at this 
point on the profile. The trajectory must have already been computed at least once prior 
to calling this routine. 



Figure 1. Basic functions. 
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Delete VTCPs: Delete the VTCPs. Remove all special vertical flags. 


T 



Compute TCP latitude and longitude data : Compute the altitude and longitude data for 
the altitude, speed, and Mach / CAS TCPs. 



Figure 1 (continued). Basic functions. 
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Algorithm Input Data 

The algorithm takes as input a list of waypoints, their trajectory-specific data, and associated 
wind profile data. The list order must begin with the first waypoint on the trajectory and end with 
the runway threshold waypoint. The trajectory-specific data includes: the waypoint's name and 
latitude / longitude data, e.g., Latitude 2 and Longitude 2 \ an altitude crossing restriction, if one 
exists, and its associated crossing angle, e.g., Crossing Altitude 2 and Crossing Angles and a speed 
crossing restriction (Mach or CAS), if one exists, and its associated CAS rate, e.g., Crossing 
CAS 2 and Crossing Rate 2 . A value of 0 as an input for an altitude or speed crossing constraint 
denotes that there is no constraint at this point. A Crossing Mach may not occur after any non- 
zero Crossing CAS input. The units for Crossing Rate are knots per second. 

For the descent from the cruise altitude, a Mach value may be specified that is different from 
the cruise Mach value. A CAS value may also be specified for the Mach / CAS transition speed 
during the descent. Additionally, an CAS speed limit at a defined altitude may also be included. 

In the U.S., this would typically be set to 250 kt at 10,000 ft. 

For the wind forecast, a minimum of two altitude reports (altitude, wind speed, and wind 
direction) should be provided at each waypoint. The altitudes should span the estimated altitude 
crossing at the associated waypoint. The algorithm assumes that the input data are valid. 

Internal Algorithm Variables 

The significant variables computed by this algorithm are: 

Altitude the computed altitude at the TCP 

CAS the computed CAS at the TCP 

DTG the computed, cumulative distance from the runway 

Ground Speed the computed ground speed at the TCP 

Ground Track the computed ground track at the TCP 

Mach the computed Mach at the TCP 

TTG the computed, cumulative time from the runway 

Additionally, the algorithm denotes TCPs in accordance with how they are generated. TCPs are 
identified as: input, from the input waypoint data; turn-entry, identifying a TCP that marks the 
start of a turn; turn-exit, identifying a TCP that marks the end of a turn; vertical TCPs (VTCPs), 
denoting a change in the altitude or speed profile; and a Mach / CAS TCP, denoting the Mach / 
CAS transition point. TCPs are also denoted relative to the associated speed value, whether the 
crossing speed is Mach or CAS derived. 
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Description of Major Functions 

The functions shown in figure 1 are described in detail in this section. The functions are 
presented in the order as shown in figure 1. Secondary functions are described in a subsequent 
section. In these descriptions, the waypoints, which are from the input data and are fixed 
geographic points, are considered to be TCPs but not all TCPs are waypoints. Nesting levels in 
the description are denoted by the level of indentation of the document formatting. Additionally, 
long sections of logic may end with end of statements to enhance the legibility of the text. 


Generate Initial Tracks and Distances 

This is an initialization function that initializes the Mach Segment flag, denoting that the speed 
in this segment is based on Mach, and calculates the point-to-point distances and ground tracks 
between input waypoints. Great circle equations are used for these calculations, noting that the 
various dimensional conversions, e.g., degrees to radians, are not shown in the following text. 

Generate the initial distances, the center-to-center distances, and ground tracks between input 
waypoints 

for (i = index number of the first waypoint; i < index number of the last waypoint; i = i + 1) 
Start with setting the Mach segments flags to false. 

Mach Segmenti- false 

Compute the waypoint-center to waypoint-center distances. 

if (i = index number of the first waypoint) Center to Center Distance t = 0 

else 


Center to Center D is lancet = 

arccosine( sine(Latitudei-i) * sine(Latitude) + cosine(Latitudei-i) * 
cosine(Latitude) * cosine(Longitudei-i - Longitude ) ) 

Ground Tracks = 

arctangent2( sinefLongitudet - Longitudet-i) * cosine (Latitude), cosine(Latitudei. 
i) * sine(Latitude) - sine(Latitudei-i) * cosinefLatitude) * 
cosine(Longitudei - Longitude^ ;)) 

end of for (i = index number of the first waypoint; i < index number of the last waypoint; i = i 

+ 1 ) ' 

Now set the runway's ground track. 

Ground Lracktost waypoint Ground Lrackiost waypoint -i 
The cumulative distance, DTG, is computed as follows: 

LI TG; mt waypoint 0 
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for (i = index number of the last waypoint; i > index number of the first waypoint; i = i- 1) 
DTGj.i = DTGi + Center to Center Distancet 

Initialize Waypoint Turn Data 

The Initialize Waypoint Turn Data function is used to determine if a turn exists at a waypoint 
and if so, inserts turn-entry and turn-exit TCPs. Waypoints that have more than a 3 degree change 
in ground track between the previous waypoint and the next waypoint are considered turn- 
waypoints. The function is performed in the following manner: 

i = index number of the first waypoint + 1 

Last Track = Ground Trackf irst waypoint 

Note that the first and last waypoints cannot be turns. 

while (i < index number of the last waypoint) 

Track Angle After = Ground Track 
a = DeltaAngle(Last Track, Track Angle After) 

Check for a turn that is greater than 135 degrees. 
if (absolute(a) > 135) 

Set an error and ignore the turn. 
a = 0 

If the turn is more than 3-degrees, compute the turn data. 
if (absolute(a) > 3) 
half turn = a/ 2 

Track Angle Center = Last Track + half turn 

This is the center of the turn, e.g., the original input waypoint. 

Ground Track = Track Angle Center 
Turn Data Trackl, = Last Track 
Turn Data Track2, = Track Angle After 
Turn Data Turn Radius t = 0 
Turn Data Path Distance t = 0 
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Insert a new TCP at the end of the turn. 


The new TCP is inserted at location i+1 in the TCP list. The TCP is inserted between 
TCP, and TCP,+ / from the original list. The function InsertWaypoint should be 
appropriate for the actual data structure implementation of this function. 

InsertWaypoint(i + 1) 

Note that TCP i+ / is the new TCP. 

TCP i+ 1 = turn-exit 

DTG i+I = DTG t 

Ground Track i+ 1 = Track Angle After 
The start of the turn TCP is as follows, 

InsertWaypoint(i) 

TCPi = turn-entry 

Note that the original TCP is now at index i + 1 . 

DTGi = DTG i+ , 

Ground Track , = Last Track 
Last Track = Track Angle After 
i= i+ 2 

end of if (absolute(a) > 3) 
else Last Track = Ground Track t 
i = i + 1 

end of while (i < index number of the last waypoint) 

Effectively, this function: 

- Marks each turn-waypoint and sets its ground track angle to the computed angle at the 
midpoint of the turn. 

- Inserts a co-distance turn-entry TCP before this turn-waypoint with the ground track angle 
for this turn-entry TCP set to the value of the inbound ground track angle. 

- Inserts a co-distance turn-exit TCP after this tum-waypoint with the ground track angle for 
this turn-exit TCP set to the value of the outbound ground track angle. 
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An example illustrating the inserted turn-start and turn-end TCPs is shown in figure 2. 


Turn waypoint, Ground Track} = 105 0 

• — Turn-exit, Ground Track j = 

Ground Tracks — 90° 

Turn-entry, Ground Track }_j = 90° 


DTGp, = DTG i+ j = DTGj 


Compute TCP Altitudes 

Beginning with the last waypoint, the Compute TCP Altitudes function computes the altitudes 
at each previous TCP and inserts any additional altitude TCPs that may be required to denote a 
change in the altitude profile. The function uses the current altitude constraint (TCP, in fig. 3), 
searches backward for the previous constraint ( TCP t . 3 in fig. 3), and then computes the distance 
required to meet this previous constraint. The altitudes for all of the TCPs within this distance are 
computed and added to the data for the TCPs. If the along-path distance to meet the previous 
constraint is not at a TCP, a new altitude VTCP is inserted at this distance. An example of this is 
shown in figure 4. This function is performed in the following steps: 


Crossing Altitude^ = 11000ft 



TCP i-3 tcp i-2 tcp i- 1 TCP i 

Figure 3. Input altitude crossing constraints. 
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Crossing Altitude^ = 11000ft 


Altitude^ ~ ^560 ft 


. _ _ jlJtitude^j - 7940ft crossing Altitude l = 5300 ft 

__ Crossing Angle t =2.3° 


TCP , 


z-3 


Afew top 


TCP , 


i-2 


TCP , 


z-y 


TCP, 


Figure 4. Computed altitude profile with TCP added. 


Set the current constraint index number, cc, equal to the index number of the last waypoint, 

cc = index number of the last waypoint 

Set the altitude of this waypoint to its crossing altitude, 

Altitude cc = Crossing Altitude cc 

While (cc > index number of the first waypoint) 

Determine if the previous constraint cannot be met. 

If (Altitude cc > Crossing Altitude C( ) 

The constraint has not been made. 

If this is the last pass through the algorithm, set an error condition 
Altitude cc - Crossing Altitude cc 

Find the prior waypoint index number pc that has an altitude constraint, e.g., a crossing 
altitude (Crossing Altitude pc ft 0). This may not always be the previous (i.e., cc - 1) 
waypoint. 

Initial condition is the previous TCP. 
pc= cc - 1 

while ( (pc> index number of the first waypoint) and ( (TCP pc ft input waypoint) or 
(Crossing Altitude pc = 0)) ) pc = pc - 1 

Save the previous crossing altitude, 

Prior Altitude = Crossing Altitude pc 
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Save the current crossing altitude {Test Altitude) at TCP cc and the descent angle (Test 
Angle ) noting that the first and last waypoints always have altitude constraints and except 
for the first waypoint, all constrained altitude points must have descent angles. 

Test Altitude = Crossing Altitude cc 

Test Angle = Crossing Angle cc 

Compute all of the TCP altitudes between the current TCP and the previous crossing 
waypoint. 

k= cc 

while k> pc 

If the previous altitude has already been reached, set the remaining TCP altitudes to 
the previous altitude. 

if (Prior Altitude < Test Altitude) 

for (k= k- 1; k> pc; k = k- 1) Altitude k = Test Altitude 

Set the altitude at the last test point. 

Altitude pc - Test Altitude 

else 

Compute the distance from TCP k to the Prior Altitude using the altitude 
difference between the Test Altitude and the Prior Altitude with the Test Angle. If 
there is no point at this distance, add a TCP at that distance. 

Compute the distance dx to make the altitude. 

dx = (Prior Altitude - Test Altitude) / (6076 * tangent(Test Angle)) 

Compute the altitude z at the previous TCP. 

z = ((DTGk-i - DTGt) * 6076) * tangent (Test Angle) + Test Altitude 

If there is a TCP prior to this distance or if z is very close to the Prior Altitude, 
compute and insert its altitude. 

if ( (DTGk-i < (DTGk + dx)) or ( absolute(z - Prior Altitude) < some small value) 

) 

if (absolute(z - Prior Altitude) < some small value) Altitude k -i = Prior 
Altitude 

else Altitude k -i = z 

Check to see if the constraint has been reached, if not, set an error condition. 
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if((k-l)=pc) 


if ( absolute( Altitude pc - Crossing Altitude pc ) > 100ft ) set an error here 
Always set the crossing exactly to the crossing value. 

Altitude pc = Crossing Altitude pc 
Update the Test Altitude. 

Test Altitude = Altitude k -i 

Decrement the counter to set it to the prior TCP. 

k= k-1 

end of if ( (DTGk-i < (DTGk + dx) ) or ( absolute(z - Prior Altitude) < some small 
value) ) 

else 


The altitude constraint is reached prior to the TCP, a new VTCP will need to 
be inserted at that point. The distance to the new TCP is, 

d = DTGk + dx 

Compute the ground track at distance d along the trajectory and save it as 
Saved Ground Track. 

Saved Ground Track = GetTrajGndTrk(d) 

Insert a new VTCP at location k in the TCP list. The VTCP is inserted 
between TCP*./ and TCP* from the original list. The function InsertWaypoint 
should be appropriate for the actual data structure implementation of this 
function. 

Insert Waypoint(k) 

Update the data for the new VTCP which is now TCPk. 

DTG k = d 

Altitude k = Prior Altitude 

Add the ground track data which must be computed if the new VTCP occurs 
within a turn. The functions WptlnTurn and ComputedGndTrk are described 
in subsequent sections. 

if (WptInTurn(k)) Ground Track k = ComputedGndTrk(k, d) 
else Ground Track k = Saved Ground Track 
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Compute and add the wind data at distance d along the path to the data of 
TCP k . 


GenerateWptWindProfile(d, TCPt) 

Test Altitude = Prior Altitude 

Since TCP k , has now been added prior to pc, the current constraint counter cc 
needs to be incremented by 1 to maintain its correct position in the list. 

cc= cc+ 1 

The function loops back to while k > pc. 

Now go to the next altitude change segment on the profile. 
cc = k 

The function loops back to while cc > index number of the first waypoint. 


Copy Crossing Angles 

The Copy Crossing Angles is a simple function that starts with the next to last TCP and copies 
the subsequent crossing angle if the current TCP does not have a crossing angle. E.g., 

for (i = index number of the last waypoint - 1; i > index number of the first waypoint; i= i- 

1 ) 


if ( Crossing A ngle, = 0) Crossing Angle, = Crossing Angle i+ j 


Do Descent Mach 

The Do Descent Mach function changes the descent waypoint Mach if the descent Mach, 
MachDescentMach, is different than the cruise Mach. The function also will add any required, 
additional TCPs. 

The calling program provides as input and retains the subsequent outputs for the following 
variables: Todld, TodMach, and TodMachRate. The variable Todld is the name of the top-of- 
descent waypoint and is initialized as a null string by the calling program. Since this function may 
overwrite the Mach and speed change rate for an input waypoint, these variables allow the 
function to retain the original values for Mach and speed change rate and to then reset these 
variables to their original values prior to recomputing new values. 

If the Mach value for the first waypoint is not set, i.e., the path does not start with a Mach 
segment, or there is no defined descent Mach, MachDescentMach = 0, the function terminates. 
Otherwise, 

If the previous TOD data for an input waypoint have been changed, these data are restored to 
their original values. 

fiini = false 


13 



i=0 


The last designated Mach waypoint, 

LastMachldx = 0 

The first designated CAS waypoint, 

FirstCasIdx= 0 
TodIdx= 0 

if Todld is not an empty string, then 

while ( (i <index number of the last waypoint) and (fini = false) ) 
if (Waypoint Id t = Todld) 
fini = true 

Restore the old values. 

Crossing Mach t = TodMach 
CAS Rate t = TodMachRate 
Todld = empty string 
Find the Mach and CAS waypoints. 
fini = false 
i = 0 

while ( (i <index number of the last waypoint) and (fini = false) ) 
if (Crossing Mach t > 0) LastMachldx = i 
else if (Crossing CASj > 0) 

FirstCasIdx = i 
fini = true 
i= i+ 1 

Find the TOD waypoint and Mach. 
fini = false 
i = 0 
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while ( (i <index number of the last waypoint) and (fini= false) ) 
if ( 'Altitude , ^ Altitude. first waypoint) 

Todldx = i - 1 
fini = true 

else if (Crossing Machi > 0) 

MachAtTod = Crossing Machi 
i= i+ 1 

if (7 odldx > 0), mark this TCP as the Top-of-Descent. 

Check for errors. There cannot be a programmed descent Mach if there is a downstream 
Mach restriction. 

if ( (LastMachldx > Todldx) \ \ (FirstCasIdx <= Todldx) ) mark this as an error condition 
else 

Save the Mach values for all input waypoints so that they may be reset on subsequent 
passes back to their original input values. 

if (WaypointTodidx is an input waypoint) 

copy the name ofWaypoiont To didx intoTodld 

TodMach = Crossing Mach Todrdx 

TodMachRate = Crossing Rate To didx 

if ( (WaypointTodidx is an input waypoint) && ( Crossing Rate To didx > 0)) 

CAS Rate = Crossing Rate To didx 

else CAS Rate = 0.75 kt/ sec (a default value) 

The following is added to force a subsequent speed calculation. 

Crossing Rate To didx = CAS Rate 

If the aircraft will slow during the descent, do the following: 
if (MachAtTod >= MachDescentMach) 

Overwrite the TOD Mach value. 

Crossing MachrodUx = MachDescentMach 
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else 


This is a special case where the aircraft is accelerating to the descent Mach. 

Invoke the secondary function DoTodAcceleration. This function is described in a 
subsequent section. 

Crossing Mach To didx — MachAtTod 


Compute Mach / CAS TCP 

If a Mach-to-CAS transition is required, this functions computes the Mach / CAS altitude and 
inserts a Mach / CAS TCP. This function is only performed if the input data starts with a Mach 
Crossing Speed for the first waypoint. The function determines the appropriate Mach and CAS 
values, calculates the altitude that these values are equal, and then determines the along-path 
distance where this altitude occurs on the profile. A Mach / CAS TCP is then inserted into the 
TCP list. 

Find the last Crossing Mach and the first Crossing CAS in the list. 

First CAS = 0 

i = index number of the first waypoint 

while ( (i <index number of the last waypoint) and (First CAS = 0) ) 
if (Crossing Mach t > 0) 

Last Mach = Crossing Mach, 

Last Mach Altitude = Altitudei 
else if (Crossing CAS t > 0) 

First CAS= Crossing CAS, 

CAS Rate = CAS Rate, 
i= i+ 1 

If there is a Mach / CAS transition speed input, use this value for the First CAS value. 

if (Mach CAS Transition > 0) First CAS = Mach CAS Transition 

Compute the Mach / CAS transition altitude. 

z = ( 1.0 - (((((0.2 * ((FirstCas/661.48) 20 ) + 1.0) 35 ) -1.0)/ 

(((0.2 * (LastMach 20 ) + 1.0) 35 ) - l.O)) 019026 )) / 0.00000687535 
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For an actual implementation, it would be beneficial to check for an error at this point. If z is 
greater than the altitude associated with the Last Mach TCP or if z is less than the altitude 
associated with the First CAS TCP, then an error should be noted. 

Find where z first occurs. 

i = index number of the first waypoint + 1 

finished = false 

while ( (i < index number of the last waypoint) and (finished = false)) 
if (Altitude, > z) i = i + 1 
else finished = true 
Find the distance to this altitude. 
jc = Altitude - Altitude , 
if (x<0) ratio = 0 
else ratio = (z - Altitude) / x 
d = ratio * (DTG,., - DTG) + DTG, 

Compute the ground track at distance d along the trajectory and save it as Saved Ground 
Track. 

Saved Ground Track = GetTrajGndTrk(d) 

Insert a new TCP at location i in the TCP list. The TCP is inserted between TCP,.; and TCP, 
from the original list. The function InsertWaypoint should be appropriate for the actual data 
structure implementation of this function. 

InsertWaypoint(i) 

Mark this TCP as the Mach / CAS transition TCP. 

Add the data for this new TCP. 

Crossing Mach, = Last Mach 
Crossing CAS, = First CAS 
CAS Ratei = CAS Rate 
DTG, = d 
Altitude, = z 

Crossing Angle t = Crossing Angle i+ i 
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Ground Trackj = Saved Ground Track 


Machj = Last Mach 
CASi = First CAS 

Compute and add the wind data at distance d along the path to the data of TCP t . 
GenerateWptWindProfHe(DTGj, TCP) 

Mark all TCPs from the first TCP (TCPf, rst waypoint) to TCP t .i as Mach TCPs. 


Compute Altitude / CAS Restriction TCP 

If an altitude / CAS restriction is required, the Compute Altitude / CAS Restriction TCP 
function computes the altitude / CAS restriction point and insert an altitude / CAS TCP. This is 
the (U.S.) point where the trajectory transitions through 10,000 ft and a 250 kt restriction is 
required. This function is only performed if the previously computed flag Need lOKRes friction is 
true. The function determines the along-path distance where this altitude / CAS occurs on the 
profile. A TCP is then inserted into the TCP list at this point. The restriction values are Descent 
Crossing Altitude and Descent Crossing CAS. 

Find the first TCP that is below the Descent Crossing Altitude in the list. 

i = index number of the first waypoint 

k= i 

fini = false 

while ( (i <index number of the last waypoint) and (fini = false) ) 
if (Altitudei < Descent Crossing Altitude) 
k= i 

fini = true 
i= i+ 1 

Find the last CAS restriction prior to the first waypoint below Descent Crossing Altitude. 


i = k - 1 


fini = false 
Last CAS = 0 
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while ( (i > 0) and (fini = false) ) 
if (Crossing CAS t > 0) 

Last CAS = Crossing CAS t 
fini = true 
i = i - 1 

Determine if an altitude / CAS TCP is required. If it is, add it. 
if ( (TCP k is a Mach segment) and (Last CAS > Descent Crossing CAS) ) 
i = k; 

Find the distance to this altitude. 
x = Altitude - Altitude t 
if (x<0) ratio = 0 

else ratio = (Descent Crossing Altitude - Altitude ) lx 
d = ratio * (DTG lA - DTG) + D re- 
compute the ground track at distance d along the trajectory and save it as Saved Ground 
Track. 

Saved Ground Track = GetTrajGndTrk(d) 

Insert a new TCP at location i in the TCP list. The TCP is inserted between TCP,.; and 
TCP, from the original list. The function InsertWaypoint should be appropriate for the 
actual data structure implementation of this function. 

InsertWaypoint(i) 

Mark this TCP as the altitude / CAS restriction TCP. 

Add the data for this new TCP. 

Crossing Machi = 0 

Crossing CAS, = Descent Crossing CAS 
Use a high value, arbitrary CAS rate. 

CAS Rate, = 0.75 kt/ sec 
DTGt = d 

Altitudet = Descent Crossing Altitude 
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Crossing Angle t = Crossing Angle i+1 
Set the Mach flag for TCP t to false 
Ground Track = Saved Ground Track 
Machj = 0 

CASi = Descent Crossing CAS 

Compute and add the wind data at distance d along the path to the data of TCP t . 
GenerateWptWindProfile(DTGi, TCPj) 


Test for Altitude / CAS Restriction Requirement 

The Test for Altitude / CAS Restriction Requirement function determines if the addition of an 
altitude / CAS restriction point is required. This is the (U.S.) point where the trajectory transitions 
through 10,000 ft and a 250 kt restriction is required. This function determines the value of the 
Need lOKRes friction flag. The function can only be called after an initial, preliminary trajectory 
has been generated. The restriction values are Descent Crossing Altitude and 
Descent Crossing CAS. 

Need lOKRes friction = false 

if ( (Descent Crossing Altitude > ) and (Descent Crossing CAS> 0) ) ok = true 

else ok = false 

If we don't start above 10,000ft, skip this whole routine. 
if (ok and ( Altitude first waypoint > Descent Crossing Altitude) ) 

Find the first point below Descent Crossing Altitude 
fini = false 
i = 0 

while ( (i <index number of the last waypoint) and (fini = false) ) 
if (Altitudei < Descent Crossing Altitude) 

Find the distance to this altitude, 
x = Altitudei-i - Altitudei 
if (x<0) ratio = 0 

else ratio = (Descent Crossing Altitude - Altitude) / x 
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s = ratio * ( CAS UI - CAS) + CAS t 


if (s> (Descent Crossing Cas + 2.0) ) Needl OKRestriction = true 
fini = true 
i— i+ 1 


Compute TCP Speeds 

The Compute TCP Speeds function is similar to Compute TCP Altitudes in its design. 
Beginning with the last waypoint, this function computes the Mach or CAS at each previous TCP 
and inserts any additional speed TCPs that may be required to denote a change in the speed 
profile. The function uses the current speed constraint, searches backward for the previous 
constraint, and then computes the distance required to meet this previous constraint. The speeds 
for all of the TCPs within this distance are computed and added to the data for the TCPs. If the 
along-path distance to meet the previous constraint is not at a TCP, a new speed VTCP is inserted 
at this distance. This function invokes two secondary functions, described in the subsequent text, 
with the invocation dependent on the constraint speed, whether it is a Mach or a CAS value. This 
function is performed in the following steps: 

Set the current constraint index number, cc, equal to the index number of the last waypoint, 
cc = index number of the last waypoint 

The speed of the first waypoint is set to its crossing speed. 

if (Crossing Mach first mypomt > 0) 

Mach fi rs t waypoint Crossing Machfirst waypoint 

CAS first waypoint MachToCasfMach fi^t waypoint* Altitude first waypoint) 

else 

CA S first -waypoint ^ OSS iflg CA Sfi rs t waypoint 


Mach first waypoint CasToM^QcH^CAS fi rs t waypoint* Altitude first waypoint) 

The speed of the last waypoint is set to its crossing speed, 

CAS ^ = Crossing CAS cc . 

A flag signifying that Mach segment computation has begun is set to false, 
Doing Mach = false 

While (cc > index number of the first waypoint) 

Set the Mach flag if the current TCP is the Mach / CAS transition point. 
if (TCP cc = Mach CAS Transition) Doing Mach = true 
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if (Doing Mach) ComputeTcpMach(cc) 
else ComputeTcpCas(cc) 

end of while cc > index number of the first waypoint 

Compute Secondary Speeds 

The Compute Secondary Speeds function adds the Mach values to CAS TCPs, the CAS values 
to Mach TCPs, and the ground speed values to all TCPs. This function is preformed in the 
following steps: 

Doing Mach = false 

Working backwards form the runway, compute the relevant speeds. 

for (i = index number of the last waypoint; i > index number of the first waypoint; i = i - 1) 
Set the flag if the current TCP is the Mach / CAS transition point. 
if (TCPi = Mach CAS Transition) Doing Mach = true 
if (Doing Mach) Cas, = MachToCas(Mach u Altitude ) 
elseMachj = CasToMach(Cast, Altitude ) 

Compute the ground track. 

if (i = index number of the first waypoint) track = Ground Track 
else if (WptInTurn(i) or (TCPi = turn-exit)) track = Ground Track 
else track = Ground Track-i 

Compute the ground speed. Compute the wind at this point. 

InterpolateWindWptAltitude(Wind Profile „ Altitude t, Wind Speed, Wind Direction) 

Ground Speedy = ComputeGndSpeedUsingTrack (Cast, track, Altitude it Wind Speed, 

Wind Direction) 

end of for (i = index number of the last waypoint; i > index number of the first waypoint; i = 

i-1) 


Update Turn Data 

The Update Turn Data function computes the turn data for each turn waypoint and modifies 
the associated waypoint's turn data sub-record. This function performs as follows: 

KtsToFps= 1.69 
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Nominal Bank Angle = 22 

index = index number of the first waypoint + 1 

while (index < index number of the last waypoint) 

Find the next input waypoint with a turn. 

while ( (index < index number of the last waypoint) and ( (TCP index f input waypoint) or 
(not Wptln Turn (index) ))) index = index + 1 

If there are no errors and there is a turn of more than 3-degrees, compute the turn data. 

if (index < index number of the last waypoint) 

Find the start of the turn. 

i = index - 1 

while (TCP t f turn-entry) i = i - 1 
start = i 

The following are all approximations and are based on a general, constant radius turn. 

The start of turn to the midpoint data is as follows, noting that the ground speeds for 
all points must be valid at this point. 

The overall distance d for this part of the turn is, 

d = DTG start - DTG index 

The special case with 0 distance between the points is, 

if (d<= 0) AvgGsFirstHalf= (Ground Speed start + Ground Speeding / 2 

else 

The overall average ground speed is computed as follows, noting that it is the 
sum of segment distance / overall distance * average segment ground speed. 

AvgGsFirstHalf= 0 

for (j = start; j < (index - 1); j = j + 1) 

dx = DTGj - DTGj+ / 

AvgGsFirstHalf= AvgGsFirstHalf + (dx / d) 

* (Ground Speedj + Ground Speedy ()/ 2 

Now, find the end of the turn. 
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i = index + 1 


while (TCPi f turn-exit) i = i + 1 
end= i 

Now, find the midpoint to the end of the turn. 

The overall distance for this part of the turn is, 

d = DTG index - DTG e „d 

Test for the special case, 0 distance between the points. 
if (d < 0) 

AvgGsLastHalf= (Ground Speed in dex + Ground Speed end ) / 2 

else 

Compute the overall average ground speed noting that it is the sum of 
segment distance / overall distance * average segment ground speed. 

AvgGsLastHalf= 0 

for (j = index; j < (end - l);j =j+ 1) 

dx =DTGj - DTGj+ / 

AvgGsLastHalf= AvgGsLastHalf + (dx/ d)* 

(Ground Speed j + Ground Speedy i) / 2 

end of for (j = index; j <= (end - 1); j =j + 1) 

end of else if (d < 0) 

The general equation is turn rate = c tan(bank angle) / v. If the bank angle is a 
constant, turn rate = cO / v. The Nominal Bank Angle = 22 degrees. 

cO = 57.3 * 32.2 / KtsToFps * tangent(Nominal Bank Angle) 

full turn = DeltaAngle(Ground Track starb Ground Track etu ) 

half turn = full turn / 2 

Compute the outputs from the average ground speed. 

Average Ground Speed = (AvgGsFirstHalf + AvgGsLastHalf) / 2 
Save the ground speed data in the turn data for this waypoint. 
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Turn Data Average Ground Speeding = Average Ground Speed 
w= cO / Average Ground Speed 
The time to make the turn is, 

Turn Data Turn Time index = absolute(full turn) / w 
The turn radius is, 

Turn Data Turn Radius in dex = (57.3 * KtsToFps * Average Ground Speed) 

/ (6076 * w) 

The along-path distance for the turn is, 

Turn Data Path Distance index = absolute(full turn) * Turn Data Turn Radius index / 57.3 

Save the turn data for the first half of the turn, denoted by the " 1 " in the variable 
name. 

Turn Data Casl index = CAS star t 

Turn Data Average Ground Speed l index = AvgGsFirstHalf 
Turn Data Trackl index = Ground Track start 

The Straight Distance values are the distances from the turn-entry TCP to the 
waypoint and from the waypoint to the turn-exit TCP. See the example in figure 5. 

Turn Data Straight Distance! index = Turn Data Turn Radius * 

tangent( absolute(half turn)) 



Figure 5. Turn distances for waypoint,. 


The Path Distance values are the along-the-path distances from the turn-entry TCP to 
a point one-half way along the turn and from this point to the turn-exit TCP. See the 
example in figure 5. 

Turn Data Path Distance = absolute(half turn) * Turn Data Turn Radius iK iex / 

57.3 

w= cO / AvgGsFirstHalf 
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Turn Data Turn Timel illdex = absolute(half turn) / w 


The data for the midpoint to the end of the turn, denoted by the "2" in the variable 
name, are as follows: 

Turn Data Cas2 index = CAS eild 

Turn Data Average Ground Speed2 index = AvgGsLastHalf 
Turn Data Track2 index = Ground Track end 

The distances for the second half of the turn are the same as for the first, but their 
calculations are recomputed here for clarity. 

Turn Data Straight Distance2 index = Turn Data Turn Radius index * 

tangent( absolute(half turn)) 

Turn Data Path Distance2 index = absolute(half turn) * Turn Data Turn Radius index / 
57.3 

w= cO /AvgGsLastHalf 

Turn Data Turn Time2 index = absolute(half turn) / w 
The DTG values are as follows: 

DTG sta rt = DTG index + Turn Data Path Distancel index 
DTG end = DTGindex - Turn Data Path Distance2 in de X 

Since the turn waypoints have been moved, the wind data need to be updated for the 
new locations. 

GenerateWptWindProflle(DTG sta rb TCP start ) 

Generate Wpt WindProfile(DTG end , TCP en f) 
end of if (index < index number of the last waypoint) 
index = index + 1 

end of while (index < index number of the last waypoint) 


Delete TCPs 

The Delete TCPs function deletes the altitude, speed, and Mach / CAS TCPs. The remaining 
TCPs will only consist of input waypoints, turn-entry, and turn-exit TCPS. This function also 
removes any flags that associate any remaining TCPs with a speed or altitude change, e.g., a 
waypoint marked as the 10,000 ft, 250 kt restriction. 
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Update DTG Data 

The Update DTG Data function is performed after the turn data have been updated and the 
VTCPs have been deleted. Only input, turn-entry, and turn-exit TCPs should be in the list at this 
time. 

DTG fi rs t waypoint 0 

i = index number of the last waypoint 
while (i > 0) 

Determine if there is a turn at either end and adjust accordingly. 
if (WptInTum(i)) 


DTGj.i = DTG t + Turn Data Path Distanced 

The following is the difference between going directly from the waypoint to going 
along the curved path. 

PriorDistanceOffset = Turn Data Straight Distance /, - Turn Data Path Distance f 
else PriorDistanceOffset = 0 
Find the next input waypoint. 
nn = i - 1 

while (TCP„„ f input waypoint) nn= nn- 1 
if (WptInTurn(nn)) 

The following is the difference between going directly from the waypoint to going 
along the curved path. 

DistanceOffset = Turn Data Straight Distance2„„ - TurnData.PathDistance2„„ 

The DTG to the input waypoint is then: 

DTG„„ = (Center to Center Distance { - PriorDistanceOffset - DistanceOffset) + DTG t 
The turn-exit DTG is then, 

DTG„„+ 1 = DTG m - Turn Data Path Distance2„„ 
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else 


The next waypoint is not in a turn. 

DTG nn = Center to Center Distance , - PriorDistanceOffset + DTGj 


i = nn 


end of while (i > 0) 


Check Turn Validity 

The Check Turn Validity function is performed after the turn data have been updated and the 
VTCPs have been deleted. Only input, turn-entry, and turn-exit TCPs should be in the list at this 
time. The function simple checks that there are no turns within turns. 

for (i = index number of the first waypoint; i < index number of the last waypoint; i = i + 1) 

if (DTGi < DTGt+i) mark this as an error condition 


Compute TCP Times 

Beginning at the runway (the last waypoint), work backwards and compute the TTG to each 
TCP. 

l'l'G Indl , x num h e r of [hi' last waypoint 0 

for (i = index number of the last waypoint; i > index number of the first waypoint; i= i- 1) 
Average Ground Speed = (Ground Speed,. 1 + Ground Speedi) / 2 
x = DTGj. i - DTGj 

Delta Time = 3600 * x / Average Ground Speed 
TTGj.i = TTGj + Delta Time 


Compute TCP Latitude and Longitude Data 

With the exception of the input waypoints, the Compute TCP Latitude and Longitude Data 
function computes the latitude and longitude data for all of the TCPs. 

In Turn = false 

Past Center = false 

Last Base = index number of the first waypoint 
Next Input = index number of the first waypoint 
Turn Index = index number of the first waypoint 
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Turn is Clockwise = true 


Turn Adjustment = 0 

Base Latitude = I .atitude [ us! Base 

Base Longitude = Longitude Last Base 

for (i = index number of the first waypoint; i < index number of the last waypoint; i = i + 1) 
if (TCPi == turn-entry) 

Turn Adjustment = 0 
InTurn= True; 

Find the major waypoint for this turn. 

Next Input = i + 1 

while ( (TCPffact i npu , f input waypoint) and (Next Input < index number of the last 
waypoint) ) Next Input = Next Input + 1 

Turn Index = Next Input 

Find the center of the turn. 

a = DeltaAngle(Ground Tracks Ground Track Next lnpu ) 

x = Turn Data Turn Radius Tu rn index / cosine(a) 

if (a > 0) Turn Clockwise =true 

else Turn Clockwise = false 

if (Turn Clockwise) al = Ground Track Turn index + 90 

else al = Ground Track Tur „imiex - 90.0 

Now compute the relative latitude and longitude values. The function RelativeLatLon 
is described in a subsequent section. 

RelativeLatLong(Latitude Tur „ Longitude Tur n index, a h x )> returning Center Latitude 
and Center Longitude 

end of if (TCPi = turn-entry) 

if (In Turn) 

Turn Adjustment = 0 

if (T urn Clockwise) al = Ground Track - 90 
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else al = Ground Trackj + 90 


if (TCPj = input waypoint) 

RelativeLatLong(Center Latitude, Center Longitude, al, x), returning Turn Data 
Latitudej and Turn Data Longitudej 

Compute the location for the center of the turn. 

a2 = DeltaAngle(T urn Data Track! i. Turn Data Track2) 

if (a2 > 0) b= Ground Trackj + 90 

else b = Ground Trackj - 90 

Compute the latitude and longitude from Turn Data Latitudej, Turn Data 
Longitudej, the angle b, and the distance, Turn Data Turn RadiuSj. 

RelativeLatLon(T urn Data Latitudej, Turn Data Longitudej, b, 

Turn Data Turn Radius ), returning Turn Data Center 
Latitudej and Turn Data Center Longitudej. 

end of if (TCPj = input waypoint) 

else RelativeLatLon(Center Latitude, Center Longitude, al, Turn Data Turn 
Radius Nexl Input ), returning Latitudej and Longitude, 

if (TCPj = turn-exit) 

Turn Adjustment = Turn Data Straight Distance2 Turn Index - 
Turn Data Path Distance2 Tum /mfec 

In Turn = false 
Last Base = Next Input 
Base Latitude = LatitudeLa St BaS e 
Base Longitude = Longitude Last Base 
end of if (In Turn) 
else 

if (TCPj = input waypoint) 

Turn Adjustment = 0 
Last Base = i 

Base Latitude — Latitudeia S j Base 
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Base Longitude = Longitude^ Base 


else 


RelativeLatLong(Base Latitude, Base Longitude, Ground Tracki.i, DTG Las t Base - 
DTGj + Turn Adjustment), returning Latitude ,■ and Longitudej 

end of for (i = index number of the first waypoint; i < index number of the last waypoint; i = i 

+ 1 ) ' 


Description of Secondary Functions 

The secondary functions are listed in alphabetical order. Note that standard aeronautical 
functions, such as CAS to Mach conversions, CasToMach, are not expanded in this document but 
may be found numerous references, e.g., reference 22. It may also be of interest to include 
atmospheric temperature or temperature deviation in the wind data input and calculate the 
temperature at the TCP crossing altitudes to improve the calculation of the various speed terms. 


ComputeGndSpeedU singT rack 

The ComputeGndSpeedU singTrack function computes a ground speed from track angle (versus 
heading), CAS, altitude, and wind data. 

b = DeltaAngle(track, Wind Direction) 

if (CAS <= 0) r = 0 

else r = (Wind Speed/ CasToTas ConversionfCAS, Altitude)) * sine(b) 

Limit the correction to something reasonable. 
if (absolute(r) > 0.8) r= 0.8 * r / absolute(r) 
heading = track + arcsine(r) 
a = DeltaAngle(heading, Wind Direction) 

TAS= CasToTas Conversion(CAS, Altitude) 

Ground Speed = (Wind Speed 2 + TAS 2 -2.0* Wind Speed * TAS * cosine(a) )°' s 


ComputeGndSpeedU singMach AndT rack 

The ComputeGndSpeedUsingMachAndTrack function computes a ground speed from track 
angle (versus heading), Mach, altitude, and wind data. 

CAS= MachToCas (Mach, Altitude) 

Ground Speed = ComputeGndSpeedU singTrack 
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ComputedGndT rk 

The ComputedGndTrk function computes the ground track at the along-path distance equal to 
distance., where distance must lie between TCP,.; and TCP i+I . It is assumed that the value for 
Ground Track is invalid. The function uses a linear interpolation based on DTG t . _/ and DTG i+1 , 
with the index value i input into the function and where the distance distance must lie between 
these points. 

d= DTGj.i - DTGi+i 

if (d < 0) Ground Track = Ground Tracks 
else 


a = (1.0 - (distance - DTG,+ij / d) * DeltaAngle(Ground Track sa. Ground Tracks) 
Ground Track = Ground Tracks + a 


ComputeTcpCas 

The variable cc is passed into and out of the ComputeTcpCas function. Beginning with the last 
waypoint, this function computes the CAS at each previous TCP and inserts any additional speed 
TCPs that may be required to denote a change in the speed profile. The function uses the current 
speed constraint, searches backward for the previous constraint, and then computes the distance 
required to meet this previous constraint. The speeds for all of the TCPs within this distance are 
computed and added to the data for the TCPs. If the along-path distance to meet the previous 
constraint is not at a TCP, a new speed VTCP is inserted at this distance. Because there is no 
general closed form solution to compute distances to meet the deceleration constraints, an 
iterative technique is used in this function. This function is performed in the following steps: 

While ((cc > index number of the first waypoint) and (TCP cc f Mach CAS Transition)) 

Determine if the previous constraint cannot be met. 

If (CAS cc > Crossing CASj 

If this is the last pass through the algorithm, set this as an error condition 
CAS cc - Crossing CAS cc 

Find the prior waypoint index number pc that has a CAS constraint, e.g., a crossing CAS 

(Crossing CAS pc f 0). This may not always be the previous (i.e., cc - 1) waypoint. 

Initial condition is the previous TCP. 

pc= cc - 1 

while ( (pc> index number of the first waypoint) and (TCP pc f Mach CAS Transition) 
and (Crossing CAS pc = 0) ) pc = pc - 1 

Save the previous crossing speed, 

Prior Speed = Crossing CAS pc 
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Save the current crossing speed ( Test Speed) at TCP cc and the deceleration rate (Test 
Rate) noting that the first and last waypoints always have speed constraints and except for 
the first waypoint, all constrained speed points must have deceleration rates. 

Test Speed = Crossing CAS cc 

Test Rate = Crossing Rate cc 

Compute all of the TCP speeds between the current TCP and the previous crossing 
waypoint. 

k= cc 

•while k> pc 

If the previous speed has already been reached, set the remaining TCP speeds to the 
previous speed. 

if (Prior Speed < Test Speed) 

for (k= k- 1; k> pc; k= k- 1) 

CAS k = Test Speed 

Mach k = CasToMach(CASb Altitude/) 

Set the speeds at the last test point. 

CAS pc = Test Speed 

if (Mach pc = 0) Mach pc = CasToMach(CAS pct Altitude p ) 

else 

Estimate the distance required to meet the crossing restriction using the winds at 
the current altitude. This is a first-estimation. 

Compute the time to do the deceleration. 

t = (Prior Speed - Test Speed) / Test Rate 

Compute the wind speed and direction at the current altitude. 

InterpolateWindWptAltitude(Wind Profile *, Altitude/,, Wind Speedl, Wind 
Directionl) 

The ground track at the current point is, 
if ( WptInTurn(k :)) Track = Ground Track k 
else Track = Ground Track k ./ 
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Current Ground Speed = ComputeGndSpeedUsingTrack( Test Speed, Track, 

Altitudeb Wind Speedl, Wind Directionl) 

Compute the wind speed and direction at the prior altitude. 

InterpolateWindWptAltitude(Wind Profile k .h Altitudeb Wind Speedl, Wind 
Directionl) 

The ground speed at the prior point. 

Prior Ground Speed = ComputeGndSpeedUsingTrack(Prior Speed, GndTrackbi, 

Altitudebi, Wind Speedl, Wind Directionl) 

Average Ground Speed = (Prior Ground Speed + Current Ground Speed) / 2. 

The distance estimate, dx, is Average Ground Speed * t. 

dx = Average Ground Speed * t / 3600 

Recompute the distance required to meet the speed using the previous estimate 
distance dx. 

Begin by computing the altitude, AltD, at distance dx. 

if ( Altitude k > Altitude k ,i) AltD — Altitude k 

else AltD = (6076 * dx)* tangent( Crossing Angle/) + Altitude k 

The new distance x is DTG k + dx. 

Compute the winds at AltD and distance x. 

InterpolateWindAtDistance(AltD, x, Wind Speed2, Wind Direction2) 

The track angle at this point, with GetTrajGndTrk defined in this section: 

Track2 = GetTrajGndTrk(x) 

The ground speed at altitude AltD is then, 

Prior Ground Speed = ComputeGndSpeedUsingTrack(Prior Speed, Track2, 

AltD, Wind Speed2, Wind Direction2) 

Average Ground Speed = (Prior Ground Speed + Current Ground Speed) / 2. 

dx = Average Ground Speed * t / 3600 

If there is a TCP prior to dx, compute and insert its speed. 

If the distance is very close to the waypoint, just set the speed. 
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if ( (DTGk-i < (DTGk + dx+ some small value) ) 


if (absolute(DTGk-i - DTGk - dx) < some small value) CASk-i = Prior Speed 
else 

Compute the speed at the waypoint using v 2 = v 0 2 + 2ax to get v. 

The headwind at the end point is, 

HeadWind2 = Wind Speed2 * cosine(Wind Direction2 - Ground Tracks) 
dx = DTGk-i - DTGk 

The value of CASk-i is computed using function EstimateNextCas, 
described in this section. 

CASk-i = EstimateNextCasfT sst Speed, Current Ground Speed, Prior 

Speed, Head Wind2, Altitude^ dx, Crossing 
Rate ct ) 

Determine if the constraint is met. 
if ( (k-1) = pc) 

Was the crossing speed met within 1 kt? If not, set this as an error. 

if (absolute(CAS pc - Crossing CAS pc ) > 1.0) Mark this as an error 
condition 

Always set the crossing exactly to the crossing speed. 

CAS pc = Crossing CAS pc 
Set the test speed to the computed speed. 

Test Speed = CASk-i 

Back up the index counter to the next intermediate TCP. 
k= k-1 

end of if ( (DTGk-i < (DTGk + dx + some small value) ) 
else 

The constraint occurs between this TCP and the previous TCP. A new VTCP 
needs to be added at this point. 

The along path distance d where the VTCP is to be inserted is: 
d = DTGk + dx 
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Save the ground track value at this distance. 


Saved Ground Track = GetTrajGndTrk(d) 

Insert a new VTCP at location k in the TCP list. The VTCP is inserted 
between TCP*./ and TCP* from the original list. The function InsertWaypoint 
should be appropriate for the actual data structure implementation of this 
function. 

Insert Waypoint(k) 

Update the data for the new VTCP which is now TCP k . 

DTG k = d 

The altitude at this point is computed as follows, recalling that the new 
waypoint is TCP k : 

if (Altitude^ / > Altitude k .i) Altitude k = Altitude *_; 

else Altitude k = (6076 * dx) * tangent(Crossing Angles 1 ) + Altitude 

CAS k = Prior Speed 

Add the ground track data which must be computed if the new VTCP occurs 
within a turn. The functions WptlnTurn and ComputedGndTrk are described 
in this sections. 

if (WptInTurn(k)) Ground Track k = ComputedGndTrk(k, d) 
else Ground Track k = Saved Ground Track 

Compute and add the wind data at distance d along the path to the data of 
TCP k . 


GenerateWptWindProfile(d, TCP *) 

Test Speed = Prior Speed 

Since TCP k , has now been added prior to pc, the current constraint counter cc 
needs to be incremented by 1 to maintain its correct position in the list. 

cc= cc+ 1 

end of while k > pc. 

Now go to the next altitude change segment on the profile. 
cc= k 

end of while cc > index number of the first waypoint 
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ComputeT cpMach 

The variable cc is passed into and out of the ComputeTcpMach function. This function is 
similar to ComputeTcpCas with the exception that the computed Mach rate will need to be 
recomputed with any change of altitude. Beginning with the last Mach waypoint (the Mach 
waypoint that is closest to the runway in terms of DTG), this function computes the Mach at each 
previous TCP and inserts any additional speed TCPs that may be required to denote a change in 
the speed profile. The function uses the current speed constraint, searches backward for the 
previous constraint, and then computes the distance required to meet this previous constraint. The 
speeds for all of the TCPs within this distance are computed and added to the data for the TCPs. 
If the along-path distance to meet the previous constraint is not at a TCP, a new speed VTCP is 
inserted at this distance. Because there is no general closed form solution to compute distances to 
meet the deceleration constraints, an iterative technique is used in this function. This function is 
performed in the following steps: 

While (cc > index number of the first waypoint) 

Determine if the previous constraint cannot be met. 

If (Mach cc > Crossing Mach C( ) 

If this is the last pass through the algorithm, mark this as an error condition 
Mach cc - Crossing Mach cc 

Find the prior waypoint index number pc that has a Mach constraint, e.g., a crossing 
Mach ( Crossing Mach pc f 0). This may not always be the previous (i.e., cc - 1) waypoint. 

Initial condition is the previous TCP. 

pc= cc - 1 

finished = false 

while ( (pc> index number of the first waypoint) and (TCP pc f Mach CAS Transition) 
and ( Crossing CAS pc = 0) ) pc- pc - 1 

Save the previous crossing speed, 

Prior Speed = Crossing Mach ^ 

Save the current crossing speed (Test Speed) at TCP cc and the deceleration rate ( Test 
Rate) noting that the first and last waypoints always have speed constraints and except for 
the first waypoint, all constrained speed points must have deceleration rates. 

Test Speed = Crossing Mach cc 

Test Rate = CasToMach(Altitude co Crossing Rate cc ) 

Compute all of the TCP speeds between the current TCP and the previous crossing 
waypoint. 

k= cc 
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while k> pc 

If the previous speed has already been reached, set the remaining TCP speeds to the 
previous speed. 

if (Prior Speed < Test Speed) 

for (k= k- 1; k> pc; k= k- 1) 

Mach k = Test Speed 

CAS k = MachToCas(Machk Altitude t) 

Mark TCP k as a Mach segment. 

Set the speeds at the last test point. 

Mach pc = Test Speed 

CAS pc = MachToCas(Mach pc , Altitude pc ) 

else 

Estimate the distance required to meet the crossing restriction using the winds at 
the current altitude. This is a first-estimation. 

Compute the time to do the deceleration. 

t = (Prior Speed - Test Speed) / Test Rate 

Compute the wind speed and direction at the current altitude. 

InterpolateWindWptAltitude(Wind Profit Altitudeh Wind Speedl, Wind 
Direction 1) 

The ground track at the current point is, 
if (WptInTurn(k)) Track = Ground Track k 
else Track = Ground Tracks 

Current Ground Speed = ComputeGndSpeedUsingMachAndTrack( Test Speed, 

Track, Altitudeh Wind Speedl, Wind Direction 1) 

Compute the wind speed and direction at the prior altitude. 

InterpolateWindWptAltitudefWind Profile k -i, Altitudeh Wind Speedl, Wind 

Directionl) 

The ground speed at the prior altitude and speed is, 
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Prior Ground Speed = ComputeGndSpeedUsingMachAndTrack(Prior Speed, 

GndTrackk-i, Altitude^, Wind Speedl, Wind 
Directionl) 

Average Ground Speed = (Prior Ground Speed + Current Ground Speed) / 2. 
The distance estimate, dx, is Average Ground Speed * t. 
dx = Average Ground Speed * t / 3600 

Compute the distance required to meet the speed using the previous estimate 
distance dx. 

Begin by computing the altitude, AltD, at distance dx. 

if (Altitude k >= Altitudek-i) AltD = Altitudek 

else AltD = (6076* dx) * tangent ( Crossing Angled + Altitude k 

Compute the average Mach rate. 

MRatel = CasToMach(Crossing Rate co Altitude t) 

MRate2 = CasToMach(Crossing Rate co AltD) 

Test Rate = (MRatel + MRate2) / 2 
t = (Prior Speed - Test Speed) / Test Rate 
The new distance x is DTG k + dx. 

Compute the winds at AltD and distance x. 

InterpolateWindAtDistance(AltD, x, Wind Speed2, Wind Direction2) 

The track angle at this point, with GetTrajGndTrk defined in this section, is: 
Track2 = GetTrajGndTrk(x) 

The ground speed at altitude AltD is then, 

Prior Ground Speed = ComputeGndSpeedUsingMachAndTrack(Prior Speed, 

Track2, AltD, Wind Speed2, Wind Direction2) 

Average Ground Speed = (Prior Ground Speed + Current Ground Speed) / 2. 

dx = Average Ground Speed * t / 3600 

If there is a TCP prior to dx, compute and insert its speed. 

If the distance is very close to the waypoint, just set the speed. 
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if ( (DTGk-i < (DTGk + dx+ some small value) ) 

if (absolute(DTGk-i - DTGk - dx)< some small value) 

Machk-i = Prior Speed 
Mark TCPk as a Mach segment. 

else 

Compute the speed at the waypoint using v 2 = v 0 2 + 2ax to get v. 

The headwind at the end point is, 

HeadWind2 = Wind Speed2 * cosine(Wind Direction2 - Ground Trackk-i) 

dx = DTGk-i - DTGk 

Compute the average Mach rate. 

MRatel = CasToMach(Crossing Rate co Altitude t) 

MRate2 = CasToMach(Crossing Rate co Altitude k -i) 

Test Rate - (MRatel + MRate2) / 2 

The value of Machk-i is computed using function EstimateNextmach, 
described in this section. 

Machk-i = EstimateNextMachfTest Speed, Current Ground Speed, Prior 

Speed, Head Wind2, Altitudeh dx, Test 
Rate) 

Determine if the constraint is met. 

if ( (k-1) = pc) 

Was the crossing speed met within 0.002 Mach? If not, set this as an 
error. 

if (absolute(Mach pc - Crossing Mach pc ) > 0.002 ) 

Mark this as an error condition 

Always set the crossing exactly to the crossing speed. 

Mach pc = Crossing Mach pc 

Set the test speed to the computed speed. 

Test Speed = Machk-i 
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Back up the index counter to the next intermediate TCP. 
k= k- 1 

end of if ( (DTGk-i < (DTGk + dx + some small value) ) 
else 


The constraint occurs between this TCP and the previous TCP. A new VTCP 
needs to be added at this point. 

The along path distance d where the VTCP is to be inserted is: 
d = DTGk + dx 

Save the ground track value at this distance. 

Saved Ground Track = GetTrajGndTrk(d) 

Insert a new VTCP at location k in the TCP list. The VTCP is inserted 
between TCP*./ and TCP* from the original list. The function InsertWaypoint 
should be appropriate for the actual data structure implementation of this 
function. 

Insert Waypoint (k) 

Update the data for the new VTCP which is now TCP k . 

DTG k - d 

The altitude at this point is computed as follows, recalling that the new 
waypoint is TCP k : 

if ( Altitude *+ / > Altitude^) Altitude k = Altitude k .i 

else Altitude k = (6076 * dx)* tangent ( Crossing A ngle^ ) + Altitude 

Mach k = Prior Speed 

Mark TCP k as a Mach segment. 

Add the ground track data which must be computed if the new VTCP occurs 
within a turn. The functions WptlnTurn and ComputedGndTrk are described 
in this sections. 

if (WptInTurn(k)) Ground Track k = ComputedGndTrk(k, d) 
else Ground Track k = Saved Ground Track 

Compute and add the wind data at distance d along the path to the data of 
TCP k . 
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GenerateWptWindProfile(d, TCP 
Test Speed = Prior Speed 

Since TCPk, has now been added prior to pc, the current constraint counter cc 
needs to be incremented by 1 to maintain its correct position in the list. 

cc= cc+ 1 

end of while k > pc. 

Now go to the next altitude change segment on the profile. 
cc = k 

end of while cc > index number of the first waypoint. 

DoT od Acceleration 

The DoTodAcceleration function handles the special case when there is an acceleration to the 
descent Mach at the top-of-descent. This function is invoked from Do Descent Mach, which 
passes in the index number for the TOD waypoint, Todldx, and the Mach value at the TOD, 
MachAtTod. 

Perform an initialization of flags and counters. 

fini = false 

skip = true 

cc= Todldx 

k= cc 

Make an initial guess of the distance to the new Mach value. 

Descent Speed = MachDescentMach 

Mach Ratei = CasToMach( 0.75 kt/ sec, Altitude cc ) 

Compute the time required to do the deceleration. 
t = ( Descent Speed - MachAtTod ) / Mach Ratei 
Compute the wind speed and direction at the current altitude. 

InterpolateWindWptAltitude( Wind Profile CC) Altitudecc Wind Speed, Wind Direction ) 

Get the ground track at the current point. 

if ( WptInTurn( Waypointcc) ) track = Ground Track cc 
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else track = Ground Track cc+ , 

TOD Ground Speed = ComputeGndSpeedUsingMachAndTrack( MachAtTod, track, 

Altitude cc , Wind Speed, Wind Direction ) 

Descent Ground Speed = ComputeGndSpeed,UsingMachAndTrack( Descent Speed, track, 
Altitude co Wind Speed, Wind Direction ) 

The average ground speed is as follows: 

Average Ground Speed = ( TOD Ground Speed + Descent Ground Speed )/ 2 
The distance estimate, dx, is Average Ground Speed * t with a conversion to nm. 
dx = Average Ground Speed * t / 3600 

AltitudeD = Altitude cc - (dx* 6076 ) * tan( Altitude Crossing Angle cc ) 

Now compute better estimates, doing this twice to refine the estimation. 
for ( i = 1; i <= 2; i = i + 1 ) 
skip = false 

Determine if this distance is beyond the next downstream waypoint. 
k= cc 

d=.DTG cc - dx 

while ( ( k < ( index number of the last waypoint - 1 ))&& ( DTG *+ j> d) ) 
if ( (k != cc ) && ( Crossing Ratek > 0) ) skip = True; 
k= k+ 1 

if ( WptInTurn( Waypointk) ) track = Ground Trackk 
else track = Ground Tracks i 

Descent Ground Speed = ComputeGndSpeedUsingMachAndTrack( Descent Speed, track, 

AltitudeD, Wind Speed, Wind Direction ) 

Average Ground Speed = ( TOD Ground Speed + Descent Ground Speed ) / 2 

Mach Rate 2 = CasToMach(Crossing Ratek, AltitudeD) 

Test Rate = (Mach Rate] + Mach Rate 2 ) / 2 

Compute the wind speed and direction at the new altitude. 

InterpolateWindWptAltitude( Waypointk Altitudek Wind Speed, Wind Direction ) 
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The ground speed at the this point is: 

Descent Ground Speed = ComputeGndSpeedUsingMachAndTrack( Descent Speed, 

Ground Tracks Altitudeh Wind Speed, Wind Direction ) 

The average ground speed is: 

Average Ground Speed = ( TOD Ground Speed + Descent Ground Speed ) / 2 
The distance, dx, is: 
dx = Average Ground Speed * t / 3600 
If there is a valid deceleration point, add it. 
if (’.skip) 

Add the VTCP for the end of the TOD acceleration. 
d = DTG CC - dx 

This needs to be here so the return data are valid. 

Old Ground Track = GetTrajGndTrk( d ) 

Temporarily save the wind data at the distance d into a temporary record, SavedWind. 
GenerateWptWindProfile(d, SavedWind) 

This needs to be inserted downstream from the TOD. 

k= k+ 1 

InsertWaypoint(k) 

Mark waypoint k as a VTCP 

If waypointk is not marked as any specific type of VTCP, mark it as the TOD acceleration. 
DTG k = d 

Altitude k = Altitude cc - (dx * 6076) * tan(Crossing Angles i) 

Mach k = Descent Speed\ 

Crossing Mach k = Descent Speed 
Mark waypoint k as a Mach segment. 

Crossing Rate k = 0.75 kt/ sec 
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if (WptInTurn(k)) Ground Tracks = ComputedGndTrk(k,d) 


else Ground Track k = Old Ground Track 
Add the saved wind data to this new waypoint. 
Copy SavedWind into WindData k 
else set an error condition 


DeltaAngle 

The DeltaAngle function returns angle a, the difference between Angle 1 and Angle2. The 
returned value may be negative, i.e., -180 degrees > DeltaAngle >180 degrees. 

a = Angle2 - Anglel 

Adjust "a" such that 0>a> 360 

if (a > 180) a = a - 360 


E stimateNextCas 

EstimateNextCas is an iterative function to estimate the CAS value, CAS, at the next TCP. 
Note that there is no closed-form solution for this calculation of CAS. The input variable names 
described in this function are from the calling routine. Also, the input deceleration value must be 
greater than 0, Test Rate > 0. 

CAS= Test Speed 

Set up a condition to get at least one pass. 
d= -10* dx 

size= 1.01 * (Prior Speed - Test Speed) 
count = 0 

if ( (dx > 0) and (T est Rate > 0) ) 

Iterate a solution. The counter count is used to terminate the iteration if the distance 
estimation does reach a solution within 0.001 n.mi. 
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while ( (absolute(d - dx) > 0.001 ) && ( count < 10) ) 
if (d > dx) CAS = CAS - size 
else CAS = CAS + size 
size= size/ 2 

The estimated time t to reach this speed, 
t = ( CAS - Test Speed) / Test Rate 
The new ground speed, 

Gs2 = CasToTas Conversion(guess, Altitude) - Head Wind2 
d= ((Current Ground Speed + Gs2) / 2) * (t/ 3600) 
count = count + 1 
end of the while loop 
Limit the computed CAS, if necessary. 
if (CAS > Prior Speed) CAS = Prior Speed 


EstimateNextMach 

EstimateNextMach is an iterative function to estimate the Mach value, Mach, at the next TCP. 
Note that there is no closed-form solution for this calculation. The input variable names described 
in this function are from the calling routine. Also, the input deceleration value must be greater 
than 0, Mach Rate > 0. 

Mach = Test Speed 

Set up a condition to get at least one pass. 
d= -10* dx 

size= 1.01 * (Prior Speed - Test Speed) 
count = 0 

if ( (dx > 0) and (T est Rate > 0)) 

Iterate a solution. The counter count is used to terminate the iteration if the distance 
estimation does reach a solution within 0.001 n.mi. 

while ( (absolute(d - dx) > 0.001) && (count < 10)) 

if (d > dx) Mach = Mach - size 
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else Mach = Mach + size 


size = size / 2 

The estimated time t to reach this speed, 
t = (Mach - Test Speed) / Test Rate 
The new ground speed, 

CAS= MachToCas(Mach, Altitude) 

Gs2 = CasToTas Conversion(CAS, Altitude) - Head Wind2 
d= ((Current Ground Speed + Gs2) / 2) * (t/ 3600) 
count = count + 1 
end of the while loop 
Limit the computed Mach, if necessary. 
if (Mach > Prior Speed) Mach = Prior Speed 

GenerateWptWindProfile 

The function GenerateWptWindProfile is used to compute new wind profile data. This function 
is a double-linear interpolation using the wind data from the two bounding input waypoints to 
compute the wind profile for a new VTCP, TCP*. The interpolations are between the wind 
altitudes from the input data and the ratio of the distance d at a point between TCP,./ and TCP, 
and the distance between TCP w and TCP,. E.g., 

- Find the two bounding input waypoints, TCP w and TCP,, between which d lies, e.g., 

TCPi., >d> TCPj. 

- Using the altitudes from the wind profile of TCP,, compute and temporarily save the wind 
data at these altitudes using the wind data from PCP W (e.g., Wind Speed Temporary} Aituudei) ■ 

- Compute the wind speed and wind direction for each altitude using the ratio r of the 
distances. Assuming that the difference between DTG,_/ and DTG t f 0, and that DTGi.i > 

DTGj. 

r = (DTGi.j - d) / (DTGt.i - DTG) 

Iterate the following for each altitude in the profile. 

Wind Speedy Annua,,! = ((1.0- r) * Wind Speed Tempora ry, mtudei) + (r * Wind Speedy mtudel ) 
a = DeltaAngle(Wind Direction Temp orary, Aituudei, Wind Direction Altitude) 

Wind Direction \ Aituudei = Wind Direction ^ Aiatudei + (r* a) 
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Figure 6 is an example of the computation data for a the wind computation at a 9,000 ft 
altitude. In this example, TCP UI has wind data at 10,000 and 8,000 ft and TCP, has wind data at 
9,000 ft. 


Wind Data^j, 4- 


10000ft 


Wind Datai / 
8000ft 


/ 

/ 


(1 - r) * Wind Data^j 
computed for 9000 ft 


DTG t _ i -d 


r * Wind Data ,■ 


DTG[_i - DTGj 


Wind Data; 
9000 ft ’ 


TCP i_] TCP k TCP i 

Figure 6. Example for computing a single wind data altitude. 


GetT rajectoryData 

The GetTrajectoryData function computes the trajectory data at the along-path distance equal 
to d and saves these data in a temporary TCP record. The function uses a linear interpolation 
based on the DTG values of the two TCPs bounding this distance and the distance d to compute 
the trajectory data at this point. 


GetTrajGndTrk 

The GetTrajGndTrk function computes the ground track at the along-path distance, distance, 
if ( (distance < 0) or (distance > DTGf irstwaypoinl ) ) Ground Track = Ground Trackfi rstwaypoint 

else 


Find where distance is on the path. 

i = index number of the last waypoint 

while (distance > DTGj i = i -1 

if (distance = DTG) Ground Track = Ground Track t 

else 


x — DTGi - DTGi+ / 

if(x < 0.0) r=0 

else r = (distance - DTG i+ ) / x 

dx = r* DeltaAngle(Ground Tracks Ground Track i+ ) 

Ground Track = Ground Track + dx 
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Interpolate WindAtDistance 

The function InterpolateWindAtDistan.ee is used to compute the wind speed and direction at an 
altitude, Altitude, for a specific distance, Distance, along the path. This function is a linear 
interpolation using the wind data from the input waypoints that bound the along-path distance. 

Find the bounding input waypoints. 

10 = index number of the first waypoint 

j=0; 

fini = false 

if (Distance < 0) Distance = 0 

while ( (fini = false) && (j < (index number of the last waypoint - 1) ) ) 
if ( (TCPj = input waypoint) and (DTGj >= Distance) ) iO=j 
if (DTGj < Distance ) fini = true 
end of the while loop 

11 = i0+ 1 

j = U 

fini = false 

while ( (fini = false) && (j < index number of the last waypoint) ) 
if ( (TCPj = input waypoint) and (DTGj <= Distance) ) 

H =j 

fini = true 
end of if 

J=J+1 

end of the while loop 

if (il > index number of the last waypoint) il = index number of the last waypoint 

if (iO = il) InterpolateWindWptAltitude(TCP i0 , Altitude) 

else 

Interpolate the winds at each waypoint. 

InterpolateWindWptAltitudefTCPio, Altitude), returning SpdO and DirO 
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InterpolateWindWptAltitude(TCPn, Altitude), returning Spdl and Dir 1 


Interpolate the winds between the two waypoints. 
r = (DTGw - Distance) / (DTGw - DTGu) 

Wind Speed =( (1.0 -r)* SpdO) +(r* Spdl) 
a = DeltaAngle(DirO, Dirl) 

Wind Direction = DirO + (r* a) 


Interpolate WindWptAItitude 

The function InterpolateWindWptAltitude is used to compute the wind speed and direction at 
an altitude, Altitude, for a specific TCP. This function is a linear interpolation using the wind data 
from the current TPC. 

Find the index numbers, pO and pi, for the bounding altitudes. 

p0= 0 
pi = 0 

for (k- 1; k <= Number of Wind Altitudes k = k + 1) 
if (Wind Altitude ,, k <= Altitude) p0= k 
if ( (Wind Altitude, k >= Altitude)and (pi = 0) ) pi = k 
if (pi = 0) pi = Number of Wind Altitudes ,■ 

If Altitude = Wind Altitude p0 or if Altitude = Wind Altitude pI then the wind data from that 
point is used. Otherwise, Altitude is not at an altitude on the wind profile of TCP h i.e., z = 
Wind Altitude, h then: 

r= (Altitude - Wind Altitude p o) / (Wind Altitude p i - Wind Altitude p o) 

Wind Speed =( (1 -r) * Wind Speed p o) + (r* Wind Speedy) 
a = DeltaAngle(Wind Direction p0 , Wind Direction p] ) 

Wind Direction = Wind Direction p o + (r* a) 
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RelativeLatLon 

The function RelativeLatLon computes the latitude and longitude from input values of latitude, 
BaseLat, longitude, BaseLon, angle, Angle, and range. Range. 

if (Angle = 180) Latitude = -range / 60+ BaseLat 

else Latitude = ( (Range * cos(Angle) ) / 60) + BaseLat 

if ( (BaseLat = 0) or (BaseLat = 180) ) Longitude = BaseLon 

else if (Angle = 90) Longitude = BaseLon + range/ (60 * cos(BaseLat) ) 

else if (Angle = 270) Longitude = BaseLon - Range / (60 * cos(BaseLat) ) 

else 

rl = tangent (4 5 +0.5* Latitude) 
r2 = tangent(45 + 0.5* BaseLat) 

if ( (rl = 0) or (r2 = 0) ) Longitude = 20, just some number, this is an error, 
else Longitude = BaseLon + (180 /pi * (tangent(Angle)* (log(rl) - log(r2)))) 


WptlnTurn 

The WptlnTurn function simply determines if the waypoint is between a turn-entry TCP and a 
turn-exit TCP. If this is true, then the function returns a value of true, otherwise it returns a value 
of false. 


Summary 

The algorithm described in this document takes as input a list of waypoints, their trajectory- 
specific data, and associated wind profile data. A full 4D trajectory can then be generated by the 
techniques described. A software prototype has been developed from this documentation. An 
example of the data input and the prototype-generated output is provided in the Appendix. 
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Appendix Example Data Sets 

Input Trajectory Data 

An example input trajectory data set is provided in Table A1 . 

The descent Mach is 0.82. The Mach / CAS transition speed for this example is 310 knots. 
Note that Waypoint- 18 is the runway threshold at a 50 ft crossing height. 


Table Al. Example of trajectory input data. 


Identifier 

Latitude 

Longitude 

Crossing 

Altitude 

Crossing 

Angle 

Crossing 

CAS 

Crossing 

Mach 

Crossing 

Rate 

Waypoint-01 

31.87476 

-103.244 

37000 

0 

0 

0.78 

0 

Waypoint-02 

32.48133 

-99.8635 

0 

0 

0 

0 

0 

Waypoint-03 

32.20548 

-98.9531 

0 

0 

0 

0 

0 

Waypoint-04 

32.19398 

-98.6621 

0 

0 

0 

0 

0 

Waypoint-05 

32.17042 

-98.113 

0 

0 

0 

0 

0 

Waypoint-06 

32.15959 

-97.8777 

0 

0 

0 

0 

0 

Waypoint-07 

32.34026 

-97.6623 

0 

0 

0 

0 

0 

Waypoint-08 

32.46908 

-97.5079 

0 

0 

0 

0 

0 

Waypoint-09 

32.64444 

-97.2967 

11700 

3.0 

0 

0 

0 

Waypoint- 10 

32.71448 

-97.2119 

11000 

1.1 

0 

0 

0 

Waypoint- 1 1 

32.74948 

-97.1695 

0 

0 

0 

0 

0 

Waypoint- 12 

32.97496 

-97.1783 

0 

0 

0 

0 

0 

Waypoint- 13 

33.10724 

-97.1754 

5300 

2.3 

220 

0 

0.5 

Waypoint- 14 

33.10658 

-97.0537 

4300 

1.8 

0 

0 

0 

Waypoint- 15 

33.03645 

-97.0541 

0 

0 

0 

0 

0 

Waypoint- 16 

33.00561 

-97.0542 

2400 

3.1 

170 

0 

0.5 

Waypoint- 17 

32.95953 

-97.0544 

1495 

3.0 

127 

0 

0.75 

Waypoint- 18 

32.91582 

-97.0546 

660 

3.0 

127 

0 

0.75 
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Input Wind Data 

An example wind speed data set is provided in Table A2. 


Table A2. Example of wind speed input data. 




Wind 

Wind 

Identifier 

Altitude 

Speed 

Direction 

Waypoint-01 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-02 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-03 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-04 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-05 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-06 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-07 

0 

20 

160 


10000 

50 

240 


20000 

60 

320 


40000 

70 

330 
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Table A2 (continued). Example of wind speed input data. 


Identifier 

Altitude 

Wind 

Speed 

Wind 

Direction 

Waypoint-08 

0 

20 

160 


10000 

50 

240 


20000 

60 

330 


40000 

70 

340 

Waypoint-09 

0 

20 

160 


10000 

50 

240 


20000 

60 

330 


40000 

70 

340 

Waypoint- 10 

0 

20 

160 


10000 

50 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 11 

0 

20 

160 


10000 

50 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 12 

0 

20 

160 


10000 

50 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 13 

0 

20 

160 


10000 

50 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 14 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 
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Table A2 (continued). Example of wind speed input data. 


Identifier 

Altitude 

Wind 

Speed 

Wind 

Direction 

Waypoint- 15 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 16 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 17 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 18 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 


Output Trajectory Data 

An example of the data available from this trajectory algorithm is provided in Table A3. Not 
shown, but also available, are the latitude and longitude data for each TCP. 


Table A3. Example of the trajectory output data. 


TCP type 

Identifier 

Altitude 

Mach 

CAS 

Mach 

Segment 

Ground 

Speed 

Track 

DTG 

TTG 

Input 

Waypoint-01 

37000 

0.78 

252.5 

true 

450.7 

77.1 

366.06 

3214.8 

Turn-entry 


37000 

0.78 

252.5 

true 

450.7 

77.1 

192.89 

1831.4 

Input 

Waypoint-02 

37000 

0.78 

252.5 

tme 

469.9 

93.3 

190.64 

1813.8 

Tum-exit 


37000 

0.78 

252.5 

tme 

487.5 

109.5 

188.39 

1796.9 

Turn-entry 


37000 

0.78 

252.5 

tme 

487.5 

109.5 

142.90 

1461.0 

Input 

Waypoint-03 

37000 

0.78 

252.5 

tme 

478.6 

101 

141.68 

1451.9 

Tum-exit 


37000 

0.78 

252.5 

tme 

469.1 

92.6 

140.46 

1442.6 
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Table A3 (continued). Example of the trajectory output data. 


TCP type 

Identifier 

Altitude 

Mach 

CAS 

Mach 

Segment 

Ground 

Speed 

Track 

DTG 

TTG 

Input 

Waypoint-04 

37000 

0.78 

252.5 

true 

469.1 

92.8 

126.90 

1338.6 

VTCP 


37000 

0.78 

252.5 

tme 

469.3 

93 

125.46 

1327.5 

VTCP 


36306 

0.82 

271.2 

tme 

494.5 

93 

123.28 

1311.2 

VTCP 


30337 

0.82 

310 

false 

509.6 

93 

104.53 

1176.8 

Input 

Waypoint-05 

28569 

0.793 

310 

false 

497.2 

93 

98.98 

1137.1 

Turn-entry 


25777 

0.751 

310 

false 

478.5 

93 

90.21 

1072.4 

Input 

Waypoint-06 

24818 

0.737 

310 

false 

446.6 

69.1 

87.20 

1048.9 

Tum-exit 


23858 

0.723 

310 

false 

415.4 

45.2 

84.19 

1023.8 

Input 

Waypoint-07 

19976 

0.672 

310 

false 

393.4 

45.3 

72.00 

915.2 

Input 

Waypoint-08 

16474 

0.629 

310 

false 

404.6 

45.4 

61.00 

816.0 

Input 

Waypoint-09 

11700 

0.576 

310 

false 

409.4 

45.5 

46.01 

683.4 

VTCP 


11432 

0.574 

310 

false 

408.5 

45.5 

43.71 

663.1 

Input 

Waypoint- 10 

11000 

0.524 

284.6 

false 

378.1 

45.5 

40.01 

629.3 

VTCP 


11000 

0.519 

282 

false 

375.1 

45.5 

39.65 

625.8 

Turn-entry 


10811 

0.507 

276.4 

false 

368.4 

45.5 

38.87 

618.3 

Input 

Waypoint- 11 

10382 

0.479 

262.9 

false 

340.6 

21.8 

37.12 

600.5 

VTCP 


10000 

0.453 

250 

false 

324.7 

19.3 

35.55 

583.5 

Tum-exit 


9954 

0.452 

250 

false 

308.9 

358.1 

35.36 

581.4 

Input 

Waypoint- 12 

7105 

0.429 

250 

false 

307.7 

1.1 

23.69 

445.1 

VTCP 


6474 

0.424 

250 

false 

307.3 

1.1 

21.10 

414.8 

Turn-entry 


5793 

0.391 

233.1 

false 

286.5 

1.1 

18.31 

381.0 

Input 

Waypoint- 13 

5300 

0.366 

220 

false 

270 

45.7 

16.29 

354.9 

Tum-exit 


4909 

0.363 

220 

false 

245 

90.3 

14.27 

326.6 

Turn-entry 


4556 

0.361 

220 

false 

242 

90.3 

12.42 

299.3 

Input 

Waypoint- 14 

4300 

0.359 

220 

false 

215.4 

135.3 

11.08 

278.2 

VTCP 


3987 

0.357 

220 

false 

204.1 

164.4 

10.21 

263.2 

Tum-exit 


3831 

0.35 

215.9 

false 

197 

180.3 

9.74 

254.7 

Input 

Waypoint- 15 

3009 

0.305 

191.2 

false 

170.7 

180.2 

7.24 

205.8 

Input 

Waypoint- 16 

2400 

0.268 

170 

false 

148.8 

180.2 

5.39 

164.1 

VTCP 


2140 

0.267 

170 

false 

148.9 

180.2 

4.65 

146.2 
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Table A3 (continued). Example of the trajectory output data. 







Mach 

Ground 




TCP type 

Identifier 

Altitude 

Mach 

CAS 

Segment 

Speed 

Track 

DTG 

TTG 

Input 

Waypoint- 17 

1495 

0.197 

127 

false 

105.5 

180.2 

2.62 

88.9 

Input 

Waypoint- 18 

660 

0.194 

127 

false 

106.9 

180.2 

0.00 

0.0 
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